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PREFACE 



This manual describes the 32~bit run time library (RTL) 
mathematical functions. 

Chapter 1 introduces the RTL mathematical functions. Chapter 2 

gives the needed information for all subprograms. Chapter 3 
details elementary, min and max, and miscellaneous functions. 
The appendixes present information on error messages and the 
subprograms documented in the manual. Arithmetic tables are also 
included. 

The R0 3 revision of this manual provides more information on the 
exponentiation routines. 

For a survey of the Run Time Library and more information on the 
construction of programs in the Run Time Library, refer to the 
following Perkin-Elmer publication: 

• FORTRAN VII Run Time Library (RTL) Introduction and Overview 
Reference Manual, 

Publication Number 29-578. 

Additional information applicable to the Run Time Library is 
contained in the following Perkin-Elmer publications: 

• FORTRAN VII Run Time Library (RTL) Language Extensions 
Reference Manual, 

Publication Number 29-580 

• FORTRAN VII Run Time Library (RTL) Real Time Extensions 
Reference Manual, 

Publication Number 29-581 

• FORTRAN VII Reference Manual, 
Publication Number 48-017 

The following manuals provide detailed information on related 
s of tware : 

• Common Assembly Language (CAL) Programming 
Re ference Manual, 

Publication Number 29-640 
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• OS/32 Library Loader Reference Manual, 
Publication Number 48-020 

• OS/32 Operator Reference Manual, 
Publication Number 29-574 

• OS/ 3 2 Programmer Reference Manual, 
Publication Number 29-613 
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CHAPTER 1 
I tJT RODUCTIOH 

This manual describes the mathematical subprogram Perkin-Elmer 
provides in the FORTRAN VII Run Time Library. Subprograms 
described include the Intrinsic and Basic External functions 
required by ANSI Standard X3J3/77 and the arithmetic routines 
which the compiler calls. See the FORTRAN VII Run Time Library 
Language Extensions Manual, Publication Number 29-580, for the 
subprograms which perform bit and byte manipulation and logical 
operations. 



This manual 
subprog ram: 



provides the following information for each 



1. The purpose served by the subprogram. 

2. The types of arguments and values. 

3. The range of admissible arguments. 

4. The accuracy (where relevant) of the result. 

5. Error messages and possible interrupts. 

6. A brief description of the algorithm. 

The subprograms can be called by assembly language code. The 
user must be aware of which arguments are passed to subprograms 
and which values are returned. See the FORTRAN VII Run Time 
Library Introduction and Overview, Publication Number 29-578. 

Chapter 2 contains information pertinent to all of the 
subprograms. Definitions are given to fix the terminology. 
Accuracy of the different subprograms, ranges of admissible 
values, and the generation of error messages are described. 

The Run Time Library provides descriptions of the function 
subprograms. These descriptions include external functions 
called by explicit statements in FORTRAN, intrinsic functions 
which can be called explicitly or implicitly (e.g., type 
conversion), and the involution routines called as a result of 
compiler-generated code. The programs are grouped as follows: 



Elementary functions 

Square root: 
Exponentia 1 : 
Trigonometric : 

Logarithm : 

Inverse Trigonometric 

Hyperbolic : 



SQ RT , DSQRT , CS QRT , CDSQRT 

EXP,DEXP,CEXP,CDEXP 

SIN,DSIN,CSIN,CDSIN,COS,DCOS,CCOS ,CDCOS , 

TAN,DTAN 

ALOCALOGIO, DLOG, DLOGIO, CLOG, CDLOG 

AC OS ,ASIN, ATAN,ATAN2, CACOS, DAS IN, DAT AN, 

DAT AN 2 

SINK, COS H,TANH,DSINH,DCOSH,DTANH 
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Min and Max functions: 



AMAXO 


AMAXl 


AMINO 


MINO 


MAXl 


AMINl 


DMAXl 


DMINl 


MAXO 


MINI 



Type conversion subprograms: 



DBLE 
IFIX 



FLOAT 
INT 



IDINT 
I NT 2 



Miscellaneous functions: 
Remainder 

AMOD DMOD MOD 

Manipulation of complex numbers 

AIMAG DIMAG CMPLX 

Absolute value 

ABS CABS DABS 

Transfer of sign 

DSIGN ISIGN SIGN 

Positive difference 

DIM IDIM DDIM 

Truncation and rounding 

AINT DINT ANINT 



DEL OAT 
SNGL 



CONJG 



DC ABS 



DNINT 



REAL 



DCONJG 



lABS 



NINT 



IDNINT 



Double-precision product of single-precision factors 
DPROD 
Involution (raising to powers): 



.IXXI 


.IXXR 


.IXXD 


.IXXC 


.IXXCD 


. RXXI 


.RXXR 


• RXXD 


.RXXC 


.RXXCD 


. DXXI 


.DXXR 


.DXXD 


.DXXC 


.DXXCD 


• CXXI 


• CXXR 


.CXXD 


.CXXC 


.CXXCD 


. CDXXI 


.CDXXR 


.CDXXD 


.CDXXC 


.CDXCD 



Complex arithmetic: 
.CXMPY 



,CDMPY 



.CXDIV 



.CDDIV 



1-2 
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CHAPTER 2 
GENERAL CONS IDERATIOWS 

2.1 TERMS RELATED TO ACCURACY 

Several methods describe the accuracy of a computation in 
floating-point arithmetic. In this manual/ the following terms 
are used: 



Absolute Error; 



Relative Error: 



is the absolute value of the difference 
between the true value of a desired number and 
the value as computed. 



is the ratio between the 
the true value. 



absolute error and 



Generally/ the relative error is the more useful of the two and 

better adapted to floating-point representation of numbers. To 

refer to a result as accurate "to five decimal places" is 
actually saying that the relative error is l.OE-5 or better. 

A statement of the relative accuracy of a computation is based on 
the assumption that the numerical data/ such as the argument of a 
function/ is perfectly accurate. 

In practice, numerical data generally has errors. Therefore, it 
is important to know the manner in which these errors can grow in 
the course of a computation. This is discussed for each function 
or group of functions under the heading "Propagation of Errors." 

Extended discussions are available for accuracy and 
error-propagation qualities of floating point programs. The 
reader can consult: 

Fike/ C.T./ Computer Evaluation of Mathematical Functions, 
Prentice-Hall, Inc./ 1968 

Hart, J.F., et al , Computer Approximations, Wiley/ 1968 

Sterbenz/ P./ Floating Point Computation/ Prentice-Hall 

2.2 ACCURACY OF FUNCTION SUBPROGRAMS IN GENERAL 

The subprograms in the Run Time Library can be classified by 
their use of integer or floating-point arithmetic and by the 
amount of computation they involve. 

Floating-point computations: the function subprograms described 
under Elementary Functions generally return values with a 
relative accuracy better than l.OE-7 for single-precision 
functions and better than l.OE-16 for double-precision functions. 
There are exceptions to this rule, however, in the case of 
logarithms and trigonometric functions. See Paragraphs 3.2.3 and 
3.2.4 for details. 
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The floating-point involution subprograms also have a relative 
accuracy of l.OE-7 for single precision or l.OE-16 for double 
precision. The other computational subprograms are accurate to 
a few units in the last significant bit (i.e., a relative error 
of l.OE-7 or better in single precision and l.OE-16 in double 
precision). The value l.OE-7 also applies to the subroutines 
mentioned under Complex Arithmetic, Paragraph 3.7. An exception 
in the case of AMOD or DMOD is described under Remaindering, 
Paragraph 3.5.1. 

Type conversion functions are accurate to one unit in the least 
significant bit. Those subprograms using only integer arithmetic 
are absolutely accurate. The same is true of subprograms which 
essentially do no arithmetic (such as AMINl and SIGN) . 

Ranges: The Basic External and Intrinsic function subprograms in 
the Run Time Library give meaningful results for all values of 
their arguments and are subject to the following range 
restrictions : 



Square root: 
Exponential : 

Logarithm: 
Trigonometric ; 



Inverse trigonometric 



hyperbolic 



The argument must be positive or zero. 

The argument x must satisfy 
-65*alog(16) < x < 63*alog(16), i.e., 
-180.2<x<174.6 

The argument must be positive. 

The argument of SIN, COS, or TAN must 
have an absolute value no greater than 
16**6 = 16,777,216. The "argument of 
nsIN, DCOS, or DTAN must have an 
absolute value no greater than 
16**14 = 72,057,594,037,927,936. 

The argument of ASIN, ACQS, DASIN, or 
DACOS must be between -1.0 and 1.0. 

The argument of SINK, COSH, DSINH, or 
DCOSH must have an absolute value no 
greater than 175.0. 



Complex functions; 



Restrictions on the arguments of 
complex functions are described under 
the individual functions. 



Type conversion 



The argument of IDINT, INT, or IFIX 
must be between -(2**31-1) and 2**31-1 
(2**31 = 2,147,483,648). The argument 
of INT2, must be between -32,768 and 
32,767. 



■r£-2 



48-025 FOO ROO 



Misjcellaneous functions: The second argument of AMOD, DMOD , or 

MOD, must be nonzero. The functions 
DSIGN, and SIGN require the second 
argument to be non-negative when the 
first argument is zero. 



Involution routines: The integer, real, double-precision, 

and complex routines reject a zero 
base with a negative exponent. The 
real and double-precision routines 
reject a negative base. 

2.3 GEWERATION OP ERROR MESSAGES 

Even if arguments satisfy the restrictions detailed under Ranges, 
they can lead to arithmetic underflow or overflow. When these 
conditions do occur, or when it is known that they will occur, 
the Run Time Library function subprograms cause an error message 
to be printed. The same happens if inadmissible arguments are 
passed to the subprogram. Details on specific error messages are 
given for each subprogram under Description of Function 
Subprograms. The error messages are of the form: 

ERROR n IN m AT a: d 

where: n identifies the diagnostic message, 

m is the name of the subprogram, 

a is the hexadecimal address of the point in the 
user's program from which the RTL is called, and 

d is a diagnostic message. 

See Appendix 1 for operation of the error message facility and 
suggested response to error messages. 

2.4 EXTERNAL AND INTRINSIC VERSIONS OF FUNCTIONS 

Except for the involution and complex arithmetic routines, every 
mathematical function in the RTL exists in two versions, as an 
intrinsic function or as an external function. The external 
functions are identified by their purely alphabetical names. The 
arguments are passed via an argument list. The intrinsic 
functions have names preceded by an ^. They also accept 
arguments via an argument list. Some intrinsic functions have 
argument passage in registers. Their names are preceded by a 
dot. Hence, certain functions in the RTL have 3 different 
versions, e.g. COS, 9C0S and .COS. 
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CHAPTER 3 
DESCRIPTION OF FUNCTION SUBPROGRAMS 

3.1 lUTRODUCTION 

This section contains information on each function subprogram and 
it covers these topics: the purpose served by each function, the 
range of argument values which it accepts, considerations of 
accuracy and propagation of errors, an explanation of run- time 
error messages, subprograms called by a function, and a brief 
description of the algorithm is included for the more complicated 
mathematical functions. 



When a function is in both the intrinsic and the external 
versions, the external version calls the intrinsic version to 
perform the actual computation. Therefore, the references to the 
function subprograms use the names of the intrinsic versions. 
Usually, the algorithms for single precision and double precision 
are similar. To describe them concisely, the function name and 
its arguments written in lower-case letters are generic names. A 
reference for nearly all algorithms is the book by Hart, cited 
earlier m this manual. The largest number representable in 
floating-point notation is approximately 16**63 or 4.0E75. This 
number (whether m single or double precision) is referred to as 
"OMEGA". Small values of the argument are frequently a special 
case and dependent upon the level of precision. The notation 
"eps" denotes 16**(-3) in single precision and 16**(-7) in double 
precision : 

KPS=y''3E100000' in single precision 
EPS=y"'3A100000', y''00OO000O' in double precision. 

The type Integer will always mean either Integer*2 or Integer*4. 
3.2 ELEIIENTARY FUNCTIONS 



3.2.1 Square Root 

Y=SQRT (X) 
DY=DSQRT(DX) 
CWW=CSQRT (CZZ) 
CDWW=CDSQRT (CDZZ) 

Purpose : 

Range of arguments: 

Accuracy : 



Argument 

Real 

Double precision 

Complex*8 

Complex* 16 



Value 

Real 

Double precision 

Complex*8 

Complex*16 



To compute the square root of an argument 
of the appropriate type. 

For SQRT and DSQRT, the argument must 
not be negative. For CSQRT and 
CDSQRT, the argument may be any 
complex number. 

SQRT and CSQRT are accurate to seven 
decimal places or better for all 
arguments. DSQRT and CDSQRT are 
accurate to 16 decimal places or 
better for all arguments. 
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Propagation of error: The relative error in SQRT(X) is 

approximately half the relative error in 
X. The only source of instability is the 
case where the absolute error in X is 
sizable compared to X. 

Error messages: If a negative x is received, then SQRT 

and DSQRT print error messages "NEGATIVE 
ARGUMENT". The square root of the abso- 
lute value of X is returned to the called 
program. 

Subroutines called: .CSQRT calls .CABS and .SQRT 

•C DSQRT calls . CDABS and .DSQRT 

Algorithms: .DSQRT and .SQRT use the same basic 

algorithm. If x=0.0/ then no processing 
is performed. If x<0.0/ then an error 
message is printed and x is replaced by 
its absolute value. 

Let x=(16**m)*u, l/16<_u<l 
Let m=2*n+ip, ip=0 or 1 

u=v* (.25**iq) / iq=0 or 1, l/4<.v<l 

A rational approximation, wa, to sqrt(v) 
is found first (Hart, formulas 0291 or 
0293). The approximation is improved by 
Heron's formula: 

wa is replaced by (v/wa+wa)*.5 

In double precision, this formula is 
applied twice. 

Let K=2*ip-iq 

Then sqrt (x )= (16**n ) * (2**k ) *wa 

The algorithm for CSQRT and CDSQRT 
makes use of SQRT and DSQRT, 
respectively. Let sqrt (x+iy) =a+i b. 

If x=y=0, no processing is performed. 

If x=0, then a =sqrt (abs (y ) /2 ) 
b= (y/abs (y ) }a 

If y = 0, x>0, then a=sqrt(x) 

b = 
If y = 0, x<0, then a = 

b=Kqrt ( -X ) 

Let r=cabs{x+iy) 

If x>0, then a=sqrt ( (x+r ) /2) 
b=Y/(2a) 

If x<0, then b= (y/abs (y) ) sqrt ( (-x+r ) /2) 
a=y/(2b) 
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3.2.2 Exponential Functions 



Y=EXP(X) 
DY=DEXP(DX) 
CWW=CEXP(CZZ) 
CDWW=CDEXP(CDZZ) 

Purpose : 



Argument 

Real 

Double precision 

Complex*8 

Complex* 16 



Value 

Real 

Double precision 

Complex*8 

Complex* 16 



To compute exp(x) where 
floating-point number of the 
type . 



x is a 

appropriate 



Range of arguments 



Accuracy : 



Propagation of errors: 



Error conditions: 



For EXP and DEXP the argument must lie 
between -65*alog(16) and 63*alog(16). 
The interval from -180.0 to +174.0 is 
within these bounds. For CEXP and 
CDEXP, the real part must be between 
-180.0 and +174.0 and the imaginary 
part must have absolute value less 
than 16**6. 



EXP and 

accurate 

better. 

complex 

accuracy 

propagation 

Trigonometric 

3.2.3. DEXP 

accurate to 

better. 



CEXP normally give results 

to seven decimal places or 

When the imaginary part of a 

argument is very large, this 

degrades; see remarks in 

of errors under 

Functions, Paragraph 

and CDEXP give results 

16 decimal places or 



The relative error of EXP(X) is equal 
the absolute error of X. 



to 



If the argument of EXP or DEXP is greater 
than the maximum acceptable value, then a 
message to that effect is put out and the 
number OMEGA is returned. If the 
argument of EXP or DEXP is less than the 
least acceptable value, then a different 
message is put out and 0.0 is returned. 
The same messages are put out if the real 
part of the argument of CEXP or CDEXP is 
out of the acceptable range. If the 
imaginary part of the argument of CEXP or 
CDEXP is unacceptable, another message is 
put out. CEXP returns the value (0,0) if 
the real part of the argument is less 
than the lower limit or the imaginary 
part is too large. If the imaginary part 
is acceptable, but the real part is above 
the upper limit, a large complex number 
is returned. 



Subroutines called; 



.CEXP calls .EXP, .SIN, and .COS 
.CDEXP calls .DEXP, .DSIN, and .DCOS 



48-025 FOO ROO 



3-3 



Algorithm; 



1. EXP and DEXP. 



If abs(x)<eps/ then 1+x is returned. 

If abs(x)>256 then it is treated as an 

error. Otherwise, let y=x*r/ where 

r=2.0/alog(2.0) 

y=n+f/ 0<f<l 

If n>504 then x is above the upper limit; 
if n<-520 then x is below the lower 
limit. Otherwise, exp(x) can be 
evaluated without overflow as follows: 

Let g=f-0.5 
Let z=2**(g/2) 

z is found by a rational approximation 
(Hart, formulas 1080 or 1082). 

Let n=B*m+k, 0<k<7 

s=2** (0.5*k+0.25) 

s is found from a table 

Finally : 

exp(x)=(16**m) *s*z 

2. CEXP and CDEXP 

Let the two parts of the argument be 
(xr, xi) 

If -180.0<xr<174.0, let u=exp(xr) 

If xr<-180.0, let u=0 
If xr>174.0, let u=OMEGA 

If abs (xi) >16**6, let s=c=0. 
Otherwise, s=sin(xi), c=cos(xi) 

Then: 

cexp (xi, xr) = (u*c, u*s) 

3.2.3 Trigonometric Functions 

Argument Value 



Y =SIN(X) 

Y =COS(X) 
Y=TAN(X) 

DY =DSIN(DX) 

DY =DCOS(DX) 

DY =DTAN(DX) 

CWW=CSIN(CZZ) 

CWW=CCOS(CZZ) 

CDWW=CDSIN (CDZZ) 

CDWW=CDCOS ( CDZZ ) 



Real 

Real 

Real 

Double 

Doub 1 e 

Double 

Complex*8 

Complex* 8 

Complex*16 

Complex*! 6 



precision 
precision 
precision 



Real 

Real 

Real 

Double precision 

Double precision 

Double precision 

Complex*8 

Complex* 8 

Complex*! 6 

Complex*! 6 



3-4 
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Purpose 



Range of arguments : 



Accuracy 



Propagation of errors: 



Error conditions: 



To compute the trigonometric sine, 
cosine, or tangent (assuming radian 
measure) of an argument of the 
appropriate type. 

For SIN, COS, and TAN, the argument 
must have an absolute value less than 
16**6. For DSIN, DCOS , and DTAN the 
argument must have an absolute value 
less than 16**4. For CSIN, CDSIN, 
CCOS , and CDCOS, the argument must 
have a real part less in absolute 
value than 16**6 and an imaginary part 
between -175 and +175. 

Normally, SIN, COS, TAN, CSIN, and 
CCOS return values which are accurate 
to seven decimal places or better. 
DSIN, CDSIN, DCOS, and CDCOS return 
values accurate to 16 decimal places 
or better. 

The absolute error in SIN(X), COS(X), 
DSIN(X), and DCOS(X) is no greater on 
absolute value than the error in X. The 
absolute error in TAN(X) or DTAN(X) is 
the absolute error in X times 
C0S{X)**(-2) . Thus when COS(X) is small, 
TAN(X) is large and its evaluation can be 
unstable. For the complex functions, 
when the imaginary part of the argument 
differs from 0.0, the relative error in 
the value is proportioned to the absolute 
error in the argument. 

If the argument of SIN, COS, or TAN is 
greater in absolute value than 16**6, 
then an error message is put out and .0 
is returned. If the argument of DSIN, 
DCOS, or DTAN is greater in absolute 
value than 16**14, then the same error 
message is put out and 0.0 is returned. 
If the argument of TAN or DTAN is such 
that function evaluation would lead to 
division by zero, then a different error 
message is put out and OMEGA is returned. 
If the real part of the argument of CSIN 
or CCOS has absolute value greater than 
16**6, a message is put out and (0.0, 
0.0) is returned. If the real part is 
acceptable, but the imaginary part is 
greater in absolute value than 175.0, 
a message is put out and a large complex 
number is returned. 



Subroutines called: 



•CSIN and .CCOS 
•SINK, and .COSH 

.CDSIN and .CDCOS call 
.DSINH, and .DCOSH 



call .SIN, 



DSIN, 



.COS , 
DCOS, 
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Algorithms: 1. Real and double-precision functions 

If abs(x)<eps, then 

sin(x)=x/ coslx) = l, tan(x)=x 

If abs(x) is too large (see "Error 
Conditions")/ then 0.0 is returned. 

Otherwise let 

ya=x*2*c where c = 2/7r 
Set ya=k+u v/here -l<.u<l. 

Depending on k, the evaluation of the 
function of x reduces to one of these 
forms : 

pm*sin (u* Tt / A) 
pm*cos (u* 7T /A) 
tan (u* 7r/4) 
cot (u* n /A) 

where prn=l or -1 

Polynomial or fractional approximations 
are used (Hart, formulae 3040, 3603, 
3820, 3823, 4242, and 4245) to evaluate 
the above 4 formulas. If u=0 and 

cot(u* TT /4 ) is required, then an error 
:aessage is put out and OMEGA is returned. 

2. Complex functions 

Let the argument be cz=(x,y). 

If abs(x)<16**6, let 

rc=cos(x), rs=sin(x). 
Otherwise rc=rs=0 

If abs(y)<175.0, let 

hc=5cosh(y) , hs=sinh(y) ; 

If y>175.0, let 

hc=OMEGA, hs=sOMEGA; 

If y<-175.0, let 

hc=OMEGA, hs= -OMEGA 

Then : 

ccos (c z) = (rc*hc, -rs*hs) 
csin (cz) = (rs*hc, rc*hs) 
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3.2.4 Logarithmic Functions 



Y=ALOG(X) 
Y=ALOG10 (X) 
DY=DLOG{DX) 
UY=DLOG10 (DX) 
CWW=CLOG(CZZ) 
CDWW=CDLOG(CDZZ) 



Purpose : 



Range of arguments; 



Accuracy : 



Propagation of error; 



Error conditions: 



Subroutines called: 



Argumen t 

Real 

Real 

Double precision 

Double precision 

Complex*8 

Complex*16 



Value 

Real 

Real 

Double precision 

Double precision 

Complex*8 

Complex*16 



To compute the logarithm of a floating- 
point number of the appropriate type. 
ALOGIO and DLOGlO compute logarithms to 
the base 10; the others, to the base e 
(natural logarithms) . 

ALOG, ALOGIO, DLOG, and DLOGlO require 
positive (non-zero) arguments. CLOG 
and CDLOG require a non-zero complex 
argument . 

ALOG, ALOGIO, and CLOG return values with 
an absolute error of l.OE-7 or less. 
DLOG and DLOGlO return values with an 
absolute error of l.OE-16 or less. 

The absolute error of ALOG (X ) equals the 
relative error of X. The absolute error 
of ALOGIO (X) equals the relative error of 
X multiplied by ALOGIO (e ) = . 4 3 4. 

When the argument of a real or 
double-precision function is 0, an 
error message is printed and -OMEGA is 
returned. When the argument is 
negative, a different error message is 
printed and the logarithm of the 
absolute value is returned. When the 
argument of CLOG or CDLOG is (0.0, 
0.0) an error message is printed and 
(-OMEGA 0.0) is returned. 



.ALOGl branches to .ALOG, 
branches to .DLOG. 



and .DLOGl 



.CLOG calls .CABS, .ALOGl, and . ATAN2 . 



.CDLOG calls 
.DATN2. 



. CDABS, 



DLOGl, 



and 



Algorithms: 



1. ALOG, ALOGIO, DLOG, DLOGlO. 



After stage (A) 
computation for 
that for base 
passes control 
passes control to 



in the algorithm, the 
base 10 is identical to 

e. Therefore, .ALOGl 
to .ALOG, and .DLOGl 
.DLOG. 
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(A) The error conditions are checked. 
The constant c is set to alog(2) in the 
programs . ALOG and . DLOG and to alogl0(2) 
in the programs . ALOGl and .DLOGl. 

(B) Let the argument be x 

x={16**m)*f, l/16<f<l 
Let g=(2**j)*f, 0<.j<3, .5<g<,l 

If g>sqrt (. 5) , let 
n=4*m- j 
h=(g-l)/(g+l); 

If g>sqrt ( .5) , let 
n=4*m- j-1 
h=(g-.5)/(g+.5) 

The rational function of h (Hart's 
formula 2701 or 2705) is found; the 
result is the logarithm to the base 2 of 
either g or 2*g, depending on whether g 
is greater or less than sqrt(.5). If 
this number is u, then v=n+u is the 
logarithm of x to the base 2. 
Multiplication by c gives the logarithm 
of X to the desired base. 

2. CLOG and CDLOG 

The error case (0.0, 0.0) is eliminated. 
If the argument (x,y) is not zero, let 

u=alog (cabs (x, y ) ) 
v=atan2 (y,x) 



3.2.5 



Then (u,v) is returned 
Inverse Trigonometric Functions 



Y=ACOS(X) 
Y=ASIN(X) 
Y=ATAN(X) 
Y=ATAN2 (XI, X2) 
Y=DACOS (X) 
Y=DASIN (X) 
Y=DATAN (X) 
Y2=DATAN2 (XI, X2) 



ARGUMENT 



Real 
Real 
Rea i 

Real 

Double 

Double 

Double 

Double 



precision 
precision 
precision 
precision 



VALUE 

Real 

Real 

Real 

Real 

Double 

Double 

Double 

Double 



precision 
precision 
precision 
precision 
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Purpose : 



To compute the inverse sine,^ cosine, or 
tangent of a number of the given type. 
The routines ATAN2 and DATAN2 compute the 
angle between an upward vertical line and 
the point with coordinates (xl/x2). When 
x2>0, this is atan(xl/x2). In general/ 
the value of ATAN2 is between - n and TT . 



Programming note: 



The definition of the value of ATAN2 is 
that which the ANSI standard requires. 
To get a value between and 2* n t one 
can take; 



TT +atan2 ( -X 1, -x2) v/hen xl<0 



Range of arguments: 



ATAN and DATAN accept any argument of the 



correct type. ATAN2 
any pair of arguments 
0.0) . ASIN, ACQS, 
accept any arguments 
1.0 . 



and DATAN2 accept 
other than (0.0, 
DASIN, and DACOS 
between -1.0 and 



Accuracy : 



The single-precision functions are 

accurate to seven places. The 

double-precision functions are accurate 
to 16 places. 



Propagation of error: 



Both absolute and relative errors are 
decreased by ATAN and DATAN. When X is 
large, the absolute error in ATAN (X) is 
very small. For ATAN2 and DATAN2, it is 
convenient to define the relative error 
of a pair of arguments as: 



sqrt ( (erl**2+er2**2) / (xl **2+x2**2) ) 

where erl and er2 are the errors in xl 
and x2 . Then, the absolute error in 
ATAN2 or DATAN2 is no greater than the 
relative error in the arguments. The 
functions ASIN, ACQS, DASIN, and DACOS 
increase the relative error; they become 
unstable as the argument nears 1.0 or 
-1.0 . 



Error conditions: 



If ATAN2 and DATAN2 receive (0.0, 0.0) as 
arguments, they put out an error message 
and return + tt /4 . If ASIN, DASIN, ACQS, 
and DACOS receive an argument greater 
than 1.0 in absolute value, they put out 
an error message and return the value 
0.0 . 



Subroutines called: 



ASIN and .ACQS call .SQRT. 
DACOS call .DSQRT. 



.DASIN and 
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Algorithms: 1. ASIN, ACQS, DASIN, DACOS 

If abs(x)<eps, then 

asin(x)=x, acos (x) = rr /2-x 

If -.5<x<.5, then asin(x) is found by a 
rational approximation (Hart, 
formula 4693 or 4698); 
a cos (x)=7r /2-asin(x) . 

If -l<x<-.5, let 

y=sqrt (0.5*x+0.5) 
u=asin (y) 
asin(x)=- 7r/2+2*u 
a cos (x)= TT /2-2*u 

If .5<x<l, let 

y=sqrt (-0.5*x+0.5) 
u=asin (y ) 
asin (x)= TT /2-2*u 
acos (x)=2*u 

2. ATAN, DATAN 

If the argument is so large that the 
value of the function is + TT /2 to the 
limit of precision involved, then that 
value is returned. Otherwise, the range 
is reduced by the following process. The 
argument belongs to one of seven 
intervals, the dividing points being: 

-tan(75«), -tan(45°), -tan(15»), 
tandS"), tan{45°), tan(75°). The 
argument is shifted to the range 
-tan{15«')<w<tan(15° ) by setting: 

w=(x-a)/ (1+ax) 

where a=tan ( (30° )m) , with -3<m<3 (in the 
extreme cases, the transformation is 
w=-l/x) . Then atan(w) is evaluated; it 
is equal to w if abs(w)<eps, and 
otherv^rise w*P (w*w) where P is a 
polynomial. The coefficients of p are 
taken from Hart (No. 4940 for ATAN, No. 
4945 for DATAN) . Then, 

atan(x)= atan (w) +m* ( rr /6) . 
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The algorithms for ATAN2 and DATAN2 are 
like those for ATAN and DATAN with the 
exception of the following: the case 
(0.0/0.0) generates an error message as 
mentioned above; in the case (xl/ 0.0)/ 
the value +7r/2 is returned if xl>0.0 and 
- 7r/2 if xKO; in the case (0.0, x2), the 
value 0.0 is returned if x2>0 and TT if 
x2<0. 

The exponents of the floating-point 
numbers xl and/ x2 are checked; and/ if 
they differ by more than 6, the smaller 
is treated as 0.0 according to the 
previous paragraph. The exception is the 
case when xl=0.0 and x2>0; then xl/x2 is 
returned . 

In the general case/ the argument pairs 
are partitioned into sectors depending on 
the signs of the arguments and the size 
of their ratio. The evaluation is 
reduced to finding atan(w) where 
w<tan(15*). The value of atan(w) is 
computed exactly as in ATAN or DATAN. 



3.2.6 Hyperbolic Functions 



Y=SINH(X) 

Y=COSH(X) 

Y=TANH(X) 

Y=DSINH(X) 

Y-DCOSH(X) 

Y=DTANH(X) 



ARGUMENT 

Real 

Real 

Real 

Double precision 

Double precision 

Double precision 



VALUE 

Real 

Real 

Real 

Double 

Double 

Double 



precision 
precision 
precision 



Range of arguments: 



SINK, COSH/ DSINH/ and DCOSH accept 
arguments between -175.0 and +175.0. 
TANK and DTANH accept all values.. 



Accuracy ; 



The single-precision functions are 

accurate to seven decimal places; the 
double-precision functions are accurate 
to 16 decimal places. 



Propagation of error: 



SINH multiplies the absolute error in X 
by COSH(X) and COSH multiplies it by 
SINH(X). When X is greater than 1 in 
absolute value, the relative error in 
SINH(X) or COSH(X) is close to tlie 
absolute error in X. TANK decreases both 
relative and absolute errors. 



48-025 FOO ROO 



3-11 



Error conditions: 



Subroutines called; 



If the argument is greater than 175.0, 
then SINK and COSH put out an error 
message and return OMEGA. If the 
argument is less than -175.0/ then SINK 
and COSH put out an error message and 
return, respectively, -OMEGA and OMEGA. 

•SINK and .COSH call . HYEXP and .TANK 
calls .EXP. 



Algorithms: 



.DSINH and . DCOSH call .HYDEX, and . DTANH 
calls .DEXP. 



1. SINH and DSINH 

If abs(x)<eps then sinh(x)=x 

If abs(x)>175.0 then an error condition 
exists as described above. 

If -175.0<x<-.5, let 

xa = -x-a log ( 2) 
u=-exp (xa) 

If u<-16**3 in single precision or 
u<-16**7 in double precision, 

then 

s i nh ( X ) =u , 
o therv/i se 

sinh(:0 =u-. 2 5/u 

If -.5<x<.5, a rational approximation is 
used (Hart, formula 2023 or 2028). 

If .5<x<175.0 ,let 

xa=x-a log ( 2 ) 
u=exp(xa) 

If u>16**3 in single precision or 
u>16**7 in double precision, 

then 

sinh(x) =u, 
otherwise 

sinh(x) =u-. 25/u 
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2. COSH and DCOSH 

If x<0, replace x by -x 

If x>175.0, an error condition exists as 

described above. If x<eps/ 

Then 

cosh(x)=1.0 

Otherwise let 

xa=x-a log (2) 
u=exp (xa) 

If u>16**3 in single precision or u>16**7 
in double precision/ 

Then 

cosh(x) =u 
Otherwise 

cosh(x)=u+.25u 

3. TANH and DTANH 

Let top=8.5 in single precision and 
21.0 in double precision. 

If x>top/ tanh(x)=l 
If x<-top, tanh(x)=-l 
If .5*a log{2) <x<top, let 

u=exp (2*x) 

tanh(x) = (u-l)/(u+l) 

If abs(x)<eps/ then tanh(x)=x 

If eps<abs(x) <• 5*alog (2) , a rational 
approximation is used. 

If -top<x<-. 5*a log(2) , let 

u=exp {"2*x ) 
tanh{x) = (l-u)/ (1+u) 
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3.3 MIN AND MAX FUNCTIONS: 



I=MaXO (11, . 
I=MINO ( 11, . 
X=AMAXO (11, 
X=AM1N0 (11, 
I=MAX1 (XI, . 
1=MIN1 (XI, . 
X=AMAX1 (XI, 
X=AM1N1 (XI, 
DX=DMAX1 (DXl, 
DX=DM1N1 (DXl, 





ARGUMENT 


VALUE 


, In) 


Integer 


Integer*4 


,In) 


Integer 


Integer*4 


. ,ln) 


Integer 


Real 


. ,ln) 


Integer 


Real 


,Xn) 


Real 


Integer*4 


,Xn) 


Real 


Integer*4 


. ,Xn) 


Real 


Real 


. ,Xn) 


Real 


Real 


. . ,DXn) 


Double 


Double 


. . ,DXn) 


Double 


Double 



where 2<n<256 



Purpose: 



To find the algebraic minimum or maximum 
of a collection of arguments of a given 
type and (in some cases) to convert to a 
different type. Types of arguments and 
values are shown in Table 3-1. 



TABLE 3-1 TYPES OF MIN AND MAX ARGUMENTS 



VALUE 


1 

ARGUMENT 








INTEGER 


REAL 


DOUBLE 


INTEGER* 4 






PRECISION 


MAXO 


MAXl 






MINO 


MINI 




REAL 


AMAXO 


AMAXl 






AMINO 


AMINl 




DOUBLE 






DMAXl 


PRECISION 






DMINl 
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Range of arguments: Unrestricted, except for MAXl and MINI, 

which convert real numbers to integers. 
MAXl and MINI require the chosen value to 
be between -2,147,483,648 and 
+ 2,147,483,647. 

Accuracy: Absolute. 

Propagation of error: Any error in the argument chosen as 

minimum or maximum is passed through 
unchanged. 



Error messages: If MAXl or MINI, try to convert a number 

too large for the appropriate type to 
fixed point, they print ^OVERFLOW' and 
return the largest integer of the 
appropriate sign and type. 
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3.4 TYPE CONVERSIONS 



Purpose : 



To convert a number from one type to 
another. The types of arguments and 
values are shown in Table 3-2. 



Range of arguments; 



Accuracy :: 



The functions with INTEGER*2 value 
require arguments between -32,768 and 
32,767. The functions with INTEGER*4 
value require arguments between 
-2,147,483,6 48 and +2,147,483,647. 

SNGL truncates the least significant 32 
bits of the fraction in a 
double-precision number. The functions 
which convert floating point to integer 
all truncate those bits which represent 
values less than one. For large 
arguments, FLOAT also truncates some 
bits; the result is accurate to the least 
significant bit in a floating point 
number. 



Propagation of error; 



A small error in the argument can be 
entirely lost or it can result in a large 
error in the value, depending on whether 
or not it moves the argument from below 
an integer to above that integer- If the 
true value of X is 2.93, then INT(X)=2. 
If X is computed as 2.85 (an error of 
.08), then IWT(X) is correctly computed 
as 2. But if X is computed as 3.0 1 (also 
an error of .08), then IWT(X) is computed 
as 3 . 



Error messages 



If the functions with integer values 
receive an argument which is too large, 
they print "OVERFLOW" and return the 
largest integer of the appropriate sign 
and type. 



Remarks; 



The name INT2, INT, REAL, DBLE, and CMPLX 
are generic names for FORTRAN VII O and 
FORTRAN VII D. This means that INT2, 
INT, REAL, DBLE, and CMPLX can take any 
of the following types of arguments: 



- Integer*2 

- Integer*4 

- Real*4 

- Double Precision (Real*8) 

- Complex*8 

- Complex*! 6 
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Table 3-2 lists their 
V a 1 ue s . 



names, types of arguments, and return 



TABLE 3-2 TYPE CONVERSION ROUTINE NAMES, TYPES OF 
ARGUMENTS AND VALUES 









TYPE 


OF INPUT 




Generic 


Speci fie 
RTL 


FORTRAN 
User 


ARGUMENT 


TYPE 
OF 




Optional 


Name 


Name 


Name 


Arg 1 


Arg 2 


FUNCTION 




@@II2 




Integer 




In teger*2 


INT2 


@INT2 




Real 




Integer*2 




?@DI2 


INT2 


Double 




Integer*2 




@@CI2 




Complex* 8 




Integer*2 




(a@CDI2 




Complex*16 




Integer*2 


INT 


@ (3 1 INT 




Integer 




Integer 




@INT 


INT 


Real 




Integer 




@IFIX 




Real 




Integer 




@IDINT 


IFIX 


Double 




In teger 




@@CINT 


IDINT 


Complex*8 




Integer 




@@CDI 




Complex*16 




Integer 


REAL 


@REAL 


REAL 


Integer 




Real 




AFLOAT 


FLOAT 


Integer 




Real 




ia@RREA 


Real 




Real 




@SNGL 




Double 




Real 




(aiaCREA 


SNGL 


Complex*8 




Real 




@(aCDRE 




Complex*16 




Real 


DBLE 


@DFLOAT 
@DBLE 


DFLOAT 
DBLE 


Integer 
Real 




Double 
Double 




?§CDBL 


Double 




Double 




§@DDBL 




Complex*8 




Double 




@DREAL 




Complex*16 




Double 


CMPLX 


(a@ICX 




Integer 


Integer 


Complex*8 




@CMPLX 


CMPLX 


Real 


Real 


Complex*8 




@(aDCX 




Double 


Double 


Complex*8 




@@ccx 




Complex*8 


Complex*8 


Complex* 8 




@@CDCX 




Complex*16 




Complex*! 6 


DCMPLX 


@@ICD 




Integer 


Integer 


Complex*16 




@@RCD 


. 


Real*4 


Real*4 


Complex*16 




@DCMPL 


DCMPLX 


Real*8 


Real*8 


Complex*16 




@@CXCD 




Complex*8 




Complex*16 




@@CDCD 




Complex* 16 




Complex*16 



CMPLX or DCMPLX takes 2 arguments of 
the same type, if the type is not 
complex. CMPLX or DCMPLX will accept 
only 1 argument of the complex type, 
in which case CMPLX just returns that 
argument. 

A complex number CZ is stored as two real 
numbers (Zl, Z2) signifying ( Zl) + (Z2) * 
SQRT(-l). The statement CZ = COMPLX(Al, 
A2) for Al, A2 of the same type which is 
not complex, results in making Al, A2 
into 2 real numbers Zl, Z2 (by floating 
or truncation) and setting CZ = (Zl, Z2). 
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3.5 MISCELLANEOUS FUNCTIONS 
3.5.1 Remaindering 



ARGUMENT VALUE 



Y=AMOD (XI, X2) Real Real 

DY=DMOD(DXl, DX2) Double precision Double precision 

K=M0D(I1,I2) Integer Integer*4 

Purposes To find the remainder as a result of 

dividing the first argument by the 
second. 



Range of arguments: The second argument must not be zero. In 

the functions AMOD and DMOD , if the 
second argument is very small compared to 
the first, then a floating point overflow 
can occur. 

Accuracy: MOD is precise. AMOD is accurate to 

about l.OE-7 times the value of X2 and 
DMOD to about l.OE-16 times the value of 
DX2. These accuracies can be degraded 
when X1/X2 becomes large. When X1/X2 is 
greater than 16**6, the returned result 
is meaningless. 

Propagation of error: Absolute error in the first argument of 

AMOD or DMOD is passed through unchanged, 
unless it changes the integer part of 
X1/X2. In that case, it leads to an 
error of the same size as X2. Absolute 
error in the second argument is 
multiplied by the integer part of X1/X2; 
if it also causes a shift in that number, 
then an error the size of X2 is 
introduced. 

Examples : 

Let the true value of x be 10.0 and the 
true value of y, 91.0. Then a mod (y,x) = 
1.0. If y is computed as 92.1/ then 
amod (y, x) =2 . 1 ; the error is passed along. 
If y is computed as 89.9 (same error in 
opposite direction)/ then a mod (y, x ) ==9. 9; 
a large error is introduced. If y is 
computed as 91.0 and x is computed as 
10.1/ then amod (y,x)=0.1; the error ' in x 
has been multiplied by y/x. If x is 
computed as 10.2/ then amod (y,x)=9.4; a 
large error is introduced. 

Lrror messages: If the second argument is 0/ then an 

error message is printed and the value of 
the first argument is returned. 
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Manipulation of Complex Numbers 

ARGUMENT 

Y=AIMAG(CZ) Complex*8 

R=CABS(CZ) Complex*8 

CW=GONJG(CZ) Complex*8 

CDWW=DIMAG(CDZZ) Complex*16 

CDWW=CDABS(CDZZ) Complex*16 

CDWW=DCONJG(CDZZ) Complex*16 



VALUE 
Real 
Real 

Complex* 8 
Complex*16 
Complex*16 
Complex*l 6 



Purpose : 



Range of arguments: 



The complex number CZ is stored as two 
real numbers (Zl, Z2), signifying 
(Z1)+(Z2)SQRT(-1) . AIMAG and DIMAG 
return Z2. CABS and CDABS return 
SQRT( (Z1)*(Z1)+(Z2)*(Z2) ). CONJG and 
DCONJG return (Zl, -Z2) as its value. 

Unrestricted. 



Accuracy : 



Error messages: 



Not applicable except to CABS, which 
is accurate to seven decimal places; 
and CDABS to sixteen decimal places. 

CABS and CDABS can cause arithmtic 
overflow. 



Subroutines called: 



CABS calls . SQRT 
CDABS calls .DSQRT 



3.5.3 Absolute Values 



Y=ABS(X) 
[)Y = DABS (DX) 
I=IABS(J) 



ARGUMENT 

Real 

Double precision 

Integer 



VALU E 

Real 

Double precision 

Integer*4 



Purpose 



Range of arguments 



To return the absolute value of an 
appropriate argument. 

Unrestricted, except lABS return 
+2,147,483,647 for full scale negative 
number -2,147,48 3,648. 



Accuracy : 
Propagation of error: 

Error messages: 
Subroutines called: 



Wot applicable. 
Wot applicable, 

Wot applicable, 
Wone. 



48-025 FOO ROO 



3-19 



3.5.4 Transfer of Sign 



K=ISIGN( I, J) 
Z=SIGN(X,Y) 
DZ=DSIGN(DX,DY) 



ARGUMENTS 

Integer 

Real 

Double precision 



VALUE 

Integer *4 

Real 

Double precision 



Purpose : 

Range of arguments: 

Accuracy : 

Propagation of error 
Remark: 

Error messages: 
Subroutines called: 



The second argument sign is ap^/lied to 
the absolute value of the first. 

The only restriction is the first 
argument must not be if the second 
argument is negative. 

Does not apply. 

Does not apply. 

Zero is a positive number for these 
routines. 

If the first argument is and the second 
argument is negative/ then an error 
message is printed and is returned. 

None . 



Positive Difference 



K=IDIM( I, J) 
Z=DIM(X,Y) 
DZ=DDIM(DX,DY) 



ARGUMENT 

Integer 

Real 

Double precision 



VALUE 

Integer*4 

Real 

Double precision 



Purpose : 

Kange of arguments: 
Accuracy : 

Propagation of error: 
Error messages : 
Subroutines called: 



To give the value or the value of the 
first argument minus the second/ 
whichever is greater. 

Unrestricted. 

To one bit in the least significant 
place . 

Not applicable. 

Not applicable. 

None. 
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3.5.6 Truncation and Rounding 

Argument Value Function 



Truncation 



Y=AINT(X) Real Real 

Y=ANINT(X) Real Real Rounding 

DY-DINT{DX) Double precision Double precision Truncation 

DY=DNINT(DX) Double precision Double precision Rounding 

J=NINT(X) Real Integer Rounding 

J=IDNINT(DX) Double precision Integer Rounding 

To obtain integer values with truncation, use the type conversion 
routines IDINT, INT, or INT2. 

Purpose: The Truncation functions remove the 

fractional part of a number; the result 
is equal to or smaller in absolute value. 
The Rounding functions return the nearest 
integer value. A fractional part of .5 
is rounded up in absolute value; thus, 
NINT(1.5)=2 and NINT (-1 . 5) =-2 . 

^^^9^- Unrestricted for AINT/ AWINT, DINT/ and 

DNINT. 

For NINT and IDNINT, the value must be no 
greater than 2**31-1 in absolute value. 

Accuracy: Not applicable. 

Propagation of error: Not applicable. 

Error messages: An overflow condition can be incurred by 

NINT or IDNINT. 

3.5.7 Double-Precision Product of Single-Precision Factors 

Argument Value 

DY=DPR0D(X1,X2) Real Double-precision 

Purpose: The argument's product is computed as a 

double-precision number. 

Range of Arguments: That of floating-point multiplication. 

Accuracy: No error is introduced. 

Propagation of error: The relative error in the product is the 

sum of the relative errors in the 
factors . 
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Error messages; 



Multiplicative overflow or 
results in error messages. 



underflow 



Subroutines called; 



None. 



3.6 INVOLUTION (RAISING TO POWERS) 



Program Bas.e 



. IXXI 
. IXXR 
.IXXD 
. IXXC 
. IXXCD 
. RXXI 
.RXXR 
.RXXD 
. RXXC 
. RXXCD 
.DXXI 
. DXXR 
. DXXD 
.DXXC 
. DXXCD 
• CXXI 
.CXXR 
.CXXD 
. CXXC 
. CXXCD 
.CDXXI 
.CDXXR 
. CDXXD 
. CDXXC 
. CDXCD 



Integer 

Integer 

I nteger 

I n teg e r 

I nteger 

Real 

Real 

Real 

Real 

Real 

Double 

Double 

Double 

Double 

Double 

Complex 

Complex 

Complex 

Complex 

Complex 

Complex*l 6 

Complex* 1 6 

Complex*! 5 

Complex*16 

Complex*! 6 



Precision 
Prec ision 
Precision 
Precis ion 
Precision 



Exgonent 

I nteger 

Real 

Double Precision 

Complex 

Complex*16 

I nteger 

Real 

Double Precision 

Complex 

Complex*16 

Integer 

Real 

Double Precision 

Complex 

Complex*! 6 

I nteger 

Real 

Double Precision 

Complex 

Complex*! 6 

Integer 

Real 

Double Precision 

Complex 

Complex*! 6 



Value 

I nteger*4 

Real 

Double Precision 

Complex 

Complex*!6 

Real 

Real 

Double Precision 

Complex 

Complex*! 6 

Double Precision 

Double Precision 

Double Precision 

Complex*! 6 

Complex*!6 

Complex 

Complex 

Complex *! 6 

Complex 

Complex*!6 

Complex*! 6 

Complex* 16 

Complex*! 6 

Complex*! 6 

Complex*! 5 



Purpose : 



Range of arguments: 



To evaluate expressions of the form 
a**b, (for complex involutions a**b is 
defined as the principal value 

determined by EXP (b*LOG (a) ) ) . All 

arithmetic types are permissable as 
bases and exponents. The interpretation 
of exponentiation operations on mixed 
mode data types is described in the 
FORTRAN VII Reference Manual. 

These following conditions are not 
acceptable: zero base with negative 

exponent; negative base with non-integer 
exponent. 
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Accuracy! 



Propagation of error; 



.IXXl introduces no error, .RXXI, 
.RXXR, and .CXXI are accurate to six 
decimal places. .DXXI, .DXXD, and 
.CDXXI are accurate to 14 decimal 
places. 

The relative error includes: the 
relative error of the base multiplied by 
the exponent and the absolute error of 
the exponent multiplied by the natural 
logarithm of the base. 



Interfaces ; 



Arguments and values are passed by 
loading them into registers as follows: 



E.E.9SL£W^ ?35.S 



!?.?i.E°.?l§.llt Y.^.i.y.e 



Error conditions 



. I XXI 
. IXXR 
. IXXD 
. IXXC 
.IXXCD 
. RXXI 
.RXXR 
. RXXD 
. RXXC 
. RXXCD 
. DXXI 
. DXXR 
. DXXD 
.DXXC 
. DXXCD 
.CXXI 
. CXXR 
.CXXD 
. CXXC 
. CXXCD 
. CDXXI 
. CDXXR 
. CDXXD 
. CDXXC 
.CDXCD 



G13 

G13 

G13 

G13 

G13 

F14 

F14 

F14 

F14 

F14 

D14 

D14 

D14 

D14 

D14 

(F12,F14) 

(F12,F14) 

(F12,F14) 

(F12,F14) 

(F12,F14) 

(D12,D14) 

(D12,D14) 

(D12,D14) 

(D12,D14) 

(D12,D14) 



G12 
F14 
D14 

(F12,F14) 
(D12,D14) 
G13 
F12 
D14 

(F10,F12) 
(D12,D14) 
G13 
F14 
D12 

(F12,F14) 
(D10,D12) 
G13 
FIO 
D14 

(F8,F10) 
(D12,D14) 
G13 
F14 
DIO 

{F12,F14) 
(D8,D10) 



G13 
F14 
D14 

(F12,F14) 
(D12,D14) 
F14 
F14 
D14 

(F12,F14) 
(D12,D14) 
D14 
D14 
D14 

(D12,D14) 
(D12,D14) 
(F12,F14) 
(F12,F14) 
(D12,D14) 
(F12,F14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 



Unacceptable combinations of arguments 
are detected and error messages are put 
out. For a zero base and a negative 
exponent, the largest positive number is 
returned. For a negative base and 
floating-point exponent, 0.0 is 
returned. The routines .IXXI, .RXXR, 
and .DXXD can anticipate certain 

of arithmetic underflow or 
.IXXI detects all such 

internally. The routines 
IXXI can suffer exponential 



occurrences 
overflow; 
occurrences 
other than 
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underflows or overflows that cause 
arithmetic fault interrupts. In these 
cases, the user sees error messages from 
both the operating system and the Run 
Time Library. 

If either a fixed-point overflow or a 
floating-point exponential overflow is 
detected, the value returned is the 
largest number of the appropriate sign. 
I f an exponential underflow is detected, 
the value returned is 0.0. 

Algorithm: If the exponent is of integer type, the 

subroutine repeatedly squares the base 
and multiplies together those iterated 
squares that correspond to '1' bits in 
the binary expression of the exponent. 
The routines which receive real 
exponents of a single or double 
precision use the formula 
x**y=exp(y*alog(x ) ) . If the exponent is 
0, then the result 1 or 1.0 is 
invariably returned, whatever the value 
of the base. 

Remarks: These routines are called by the FORTRAN 

compiler to evaluate expressions of the 
form a**b. User can access these 
routines via assembly code. 
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APPENDIX 1 
ERROR MESSAGES 

If a mathematical RTL program detects an error/ it calls a 
routine .ERR to print out an error message in this form: 

ERROR n IN m AT a: d 

where: n is the error message number, 

m is the name of the RTL function, 

a is the hexadecimal return address in the 
user's program, and 

d is either a number or a diagnostic message. 

The messages are listed by number in the error messages list in 
this appendix. The following comments are intended to help 
interpret the messages and correct their causes. 

Messages Pertaining to the Argument List 

V'Jhen a mathematical RTL program is called by giving it an 
argument list with argument checking turned on, a series of 
messages can be put out as a result of errors in the list. 
Please refer to FORTRAN VII Run Time Library Introduction and 
Overview, Publication Number 29-578 for argument checking. 

The error messages may be one of the following: 

1. ILLEGAL CALL (72). 

A function has been called as a subroutine, e.g. 

CALL SQRT (,...). 

The solution is to revise the text of the source 
program or, if the user means to call an external 
subroutine of that name, to ensure that it is brought 
in at linkage time. 

2. INCORRECT TYPE OF VALUE OF FUNCTION (7 3). 

A function supplied in the RTL has been called as an 
external function and the user program has specified a 
value type that is incorrect. Example: 

IMPLICIT INTEGER (A-H) 
EXTERNAL COS 
A=COS(X) 

Check implicit and explicit types or include a 
user-written program of the appropriate name. 
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3. ARGUMENT OF INCORRECT CLASS (75). 

An argument has been classified as an array or 
external function and not a number. Check for the 
presence or absence of a DIMENSION statement, a 
misnamed argument/ or an incorrect EXTERNAL statement. 

4. AN ARGUMENT OF INCORRECT TYPE (45). 

The argument supplied is of the wrong type for this 
function. Check implicit and explicit type statements 
and types of expressions passed as arguments. 

5. INCORRECT NUMBER OF ARGUMENTS (76). 

Check argument lists in source program and for the 
presence or absence of commas. 

Arguments Outside Acceptable Ranges 

These messages can occur whether or not a program is called with 
arguments in registers or by giving it a list. 

1. ARGUMENT TOO LARGE (49) 
REAL PART TOO LARGE (52) 
IMAGINARY PART TOO LARGE (55) 

Tiie absolute value of the number is too large if the number is 
positive or negative. Check the numerical stability of the 
computation being performed. Check if variables have been 
initialized. Consider building error-checking logic into the 
program. 

2. ARGUMENT OUT OF RANGE: POSITIVE (5 0) 
ARGUMENT OUT OF RANGE: NEGATIVE (51) 
REAL PART OUT OF RANGE: POSITIVE (53) 
REAL PART OUT OF RANGE: NEGATIVE (54) 
IMAGINARY PART OUT OF RANGE: POSITIVE (56) 
IMAGINARY PART OUT OF RANGE: NEGATIVE (57) 

The number is greater than its maximum permissible value or less 
than its minimum permissible value. Remedies are as in #1. 

3. OVERFLOW ON COMVERSIOM (58) 

The absolute value of a floating-point number was too large to be 
converted to an integer. Besides the messages listed in 1 above, 
consider if this conversion is really appropriate. 



Al-2 48-025 FOO ROO 



APPKNDIX 1 (Continued) 

4. ZERO ARGUMENT (46) 
NEGATIVE ARGUiylENT (47) 
ARGUMENTS (0/0) (48) 
ZERO DIVISOR (62) 

ARGUMENTS (ZERO/ NEGATIVE) (6 3) 

These cases must be averted by building appropriate logic into 
the source code. Numerical instability should be closely 
checked; the difference of two numbers may be quite inaccurate. 

5. NEGATIVE EXPONENT (6 4) 
NEGATIVE BASE (74) 

ZERO BASE/ NEGATIVE EXPONENT (65) 

These are conditions that can arise in computing powers. Check 
the validity of the computations which determined the base and 
exponent; consider building in logic to avoid these erroneous 
conditions . 

6. EXPONENTIAL OVERFLOW (60) 
EXPONENTIAL UNDERFLOW (61) 

These are conditions that can arise when the computation of a 
number raised to a power has led to arithmetic faults. To 
safeguard against these conditions / build logic into the program 
to handle extreme values of either base or exponent. 
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APPENDIX 1 (Continued) 

RTL PROGRAM ERROR MESSAGES 

ERROR MESSAGES THAT MAY BE INVOKED FROM MATHEMATICAL RTL 
PROGRAMS. 

DIAGNOSTIC I^IESSAGE 

AN ARGUMENT OF INCORRECT TYPE 

ZERO ARGUMENT 

NEGATIVE ARGUMENT 

ARGUMENTS (0,0) 

ARGUMENT TOO LARGE 

ARGUMENT OUT OF RANGE : POSITIVE 

ARGUMENT OUT OF RANGE : NEGATIVE 

REAL PART TOO LARGE 

REAL PART OUT OF RANGE : POSITIVE 

REAL PART OUT OF RANGE : NEGATIVE 

IMAGINARY PART TOO LARGE 

IMAGINARY PART OUT OF RANGE : POSITIVE 

IMAGINARY PART OUT OR RANGE : NEGATIVE 

OVERFLOW ON CONVERSION 

VALUE OF SIZE IS ILLEGAL 

EXPONENTIAL OVERFLOW 

EXPONENTIAL UNDERFLOW 

ZERO DIVISOR 

ARGUMENTS (ZERO, NEGATIVE) 

NEGATIVE EXPONENT 

ZERO BASE, NEGATIVE EXPONENT 



NUMBER 


45 


46 


47 


48 


49 


50 


51 


52 


53 


54 


55 


56 


57 


58 


59 


60 


61 


62 


63 


64 


65 
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APPENDIX 1 (Continued) 

RTL PROGRAM ERROR MESSAGES (Continued) 
NUMBER DIAGNOSTIC MESSAGE 

66 NEGATIVE INDEX 

67 SIZE TOO LARGE 

68 SIZE NOT POSITIVE 

6 9 ILLEGAL FUNCTION CODE 

70 ILLEGAL TRAP CODE 

71 ILLEGAL LEVEL OB'' TRAPPING 

72 ILLEGAL CALL 

73 INCORRECT TYPE OF VALUE OF FUNCTION 

74 NEGATIVE BASE 

75 ARGUMENT OF INCORRECT CLASS 

76 INCORRECT NUMBER OF ARGUMENTS 
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APPENDIX 2 
ALPHAbtiTICAL LIST OF SUBPROGRAl'iS 
DOCUMENTED IN THIS MANUAL 



Math 


Paraqraph 


Math 


Paraqraph 


Function 


Reference 


Function 


Reference 


@ABS 


3.5.3 


@DFLOAT 


3.4 


@ACOS 


3.2.5 


@DIM 


3.5.5 


?AIMAG 


3.5.2 


(aniMAG 


3.5.2 


@AINT 


3.5.6 


(aniNT 


3.5.6 


@ALOG 


3.2.4 


@DLOG 


3.2.4 


^ALOGl 


3.2.4 


(aDLOGl 


3.2.4 


@AMAXO 


3.3 


@DMAX1 


3.3 


@AMAX1 


3.3 


laOMINl 


3.3 


@AMINO 


3.3 


@DMOD 


3.5.1 


@AMINl 


3.3 


(aONINT 


3.5.6 


@ANINT 


3.5.6 


@DPROD 


3.5.7 


?ASIN 


3.2.5 


@DREAL 


3.4 


?ATAN 


3.2.5 


@DSIGN 


3.5.4 


@ATAN2 


3.2.4 


@DSIN 


3.2.3 






@DSINH 


3.2.6 






@DSQRT 


3.2.1, 3.2.5 


(aCABS 


3.2.4 


@DTAN 


3.2.3 


(acCOS 


3.2.3 


@DTANH 


3.2.6 


@CDABS 


3.5.2 






@CDCOS 


3.2.3 


@EXP 


3.2.2., 3.2. 


@CDEXP 


3.2.2 


@ FLOAT 


3.4 


@CDLOG 


3.2.4 


@IABS 


3.5.3 


@CDSIN 


3.2.3 


@IDIM 


3.5.5 


@CDSQR 


3.2.1 


@IDINT 


3.4 


@CEXP 


3.2.2 


@IDNIN 


3.5.6 


@CLOG 


3.2.4 


@IFIX 


3.4 


@CMPLX 


3.4 


(ajNT 


3.4 


liCGNJG 


3.5.2 


(aiNT2 


3.4 


@COS 


3.2.2 


@ISIGN 


3.5.4 


@COSH 


3.2.3 






@CSIN 


3.2.3 


@MAXO 


3.3 


@CSQRT 


3.2.1 


@MAX1 


3.3 






@MINO 


3.3 






@MIN1 


3.3 


@DABS 


3.5.3 


@MOD 


3.5.1 


laDACOS 


3.2.5 


@NINT 


3.5.6 


(aOASIN 


3.2.5 






(BDATAN 


3.2.5 


@REAL 


3.4 


@DATN2 


3.2.5 


@SIGN 


3.5.4 


^DBLE 


3.4 


@SIN 


3.2.2, 3.2.3 


(aoCMPL 


3.4 


esiNH 


3 .2.3, 3.2.6 


(aDCONJG 


3.5.2 


(aSNGL 


3.4 


@DCOS 


3.2.3 


(asQRT 


3.2.1, 3.2.5 


@DCOSH 


3.2.6 






@DDIM 


3.5.5 


@TAN 


3.2.3 


@DEXP 


3.2.6 


@TANH 


3.2.6 



3.5.2 
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Math 
Function 

e@ccx 

(3(30 DBL 
@@CDCD 

@@CDI 

§?CDI2 

@@CDRE 

@@CI2 

@@CINT 

e?CREA 

ia@CXCD 

C^ 3 RC D 

• ACOS 
.AINT 
. ALOG 
. ALOG 1 
. AWINT 
.ASIN 
.ATAN 
.AT AN 2 

.CABS 
. CCOS 
. C DA BS 
.CDCOS 
.CDDIV 
. CDEXP 

.CDLOG 
. CDMPY 
. CDSIN 
. CDSQRT 
. CDXXC 
. CDXXC D 
. CDXXD 
.CDXXI 
. CDXXR 
. CEXP 
.CLOG 
.COS 
-COSH 
.CSIN 
.CSQRT 
.CXDIV 
. CXMPY 
. CXXC 
. CXXCD 
. CXXD 
. CXXR 
. CXXI 



Paraara^h 
Reference 

3.4 
3.4 
3,4 
3.4 
3.4 
3.4 
3.4 
3.4 
3.4 
3.4 
3.4 
3.4 
3.4 

3.2.5 

3.5.6 

3.2 

3.2, 

3.5, 

3. 2. 

3.2. 

3.2. 



,4 
,4 
,6 
,5 
5 
4 



3.2.4 

3.2.3 

3.2.4 

3.2. 3 

3.7 

3.2.2 



3.2 

3.7 

3. 2 

3.2, 

3.6 

3.6 

3.6 

3.6 

3.6 

3.2, 

3.2 



3.2.2, 3.2.3 
3.2.3/ 3.2.6 



3.2, 

3.2, 

3.7 

3.7 

3.6 

3.6 

3.6 

3.6 

3.7 



Math 
Function 


Paraqraph 
Reference 


. DACOS 


3.2.5 


.DASIN 


3.2.5 


. DATAN 


3.2.5 


. DATN 2 


3.2.5 


.DCOS 


3.2.3 


.DCOSH 


3.2.6 


.DEXP 


3.2.6 


.DINT 


3.5.6 


.DLOG 


3.2.4 


.DLOGl 


3.2.4 


. DNINT 


3.5.6 


.DSIN 


3.2.3 


.DSINH 


3.2,6 


. DSQRT 


3.2.1, 3 . 


.DTAN 


3.2.3 


. DTANH 


3.2.6 


.DXXC 


3.6 


.DXXCD 


3.6 


.DXXD 


3.6 


.DXXI 


3.6 


.DXXR 


3.6 



EXP 



2.5 



3.2.2, 3.2.6 



.HYDEX 


3.2.6 




.HYEXP 


3.2.6 




.IXXC 


3.6 




.IXXCD 


3.6 




.IXXI 


3.6 




.IXXD 


3.6 




.IXXR 


3.6 




.RXXC 


3.6 




. RXXCD 


3.6 




.RXXI 


3.6 




.RXXD 


3.6 




.SIN 


3.2.2, 


3.2.3 


-SINH 


3.2.3, 


3.2.6 


.SQRT 


3.2.1, 


3.2.5, 3.5.2 


.TAN 


3.2.3 




.TANH 


3.2.6 
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Math 


Paraqraph 


Function 


Reference 


ABS 


3.5.3 


ACQS 


3.2.5 


A IMAG 


3. 5.2 


A INT 


3.5.6 


ALOG 


3.2.4 


ALOGIO 


3.2.4 


AMAXO 


3.3 


AMAXl 


3.3 


AMINO 


3.3 


AMINl 


3.3 


AMOD 


3.5.1 


AN INT 


3.5.6 


A SIN 


3.2.5 


A TAN 


3.2.5 


ATAN2 


3.2.5 


CABS 


3.5.1 


CCOS 


3.2.3 


CDABS 


3.5.2 


CDCOS 


3.2.3 


CDEXP 


3.2.2 


CDLOG 


3.2.4 


CDSIN 


3.2.3 


CDSQRT 


3.2.1 


CEXP 


3.2.2 


CLOG 


3.2.4 


CME'LX 


3.4 


CONJG 


3.5. 2 


COS 


3.2.3 


COSH 


3.2.6 


CSIN 


3.2.3 


CSQRT 


3.2.1 


DABS 


3.5.3 


DACOS 


3.2.5 


DASIN 


3.2.5 


DAT AN 


3.2.5 


DATAN2 


3.2.5 


DBLE 


3.4 


DCMPLX 


3.4 


DCONJG 


3.5.2 


DCOS 


3.2.3 


DCOSH 


3.2. 6 


DDIM 


3.5.5 


DEXP 


3.2.2 


D FLOAT 


3.4 


DIM 


3.5.5 



Math 
Function 



Paraqraph 
Re ferance 



DIMAG 


3.5.2 


DINT 


3.5.6 


DLOG 


3.2.4 


DLOGIO 


3.2.4 


DMAXl 


3.3 


DMINl 


3.3 


DMOD 


3.5.1 


DNINT 


3.5.6 


DPROD 


3.5.7 


DREAL 


3.4 


DSIGN 


3.5.4 


DSIN 


3.2.3 


DSINH 


3.2.6 


DSQRT 


3.2.1 


DTAN 


3. 2.3 


DTANH 


3.2.6 



EXP 



3.2.2 



FLOAT 


3.4 


lABS 


3.5.3 


IDIM 


3.5.5 


I DINT 


3.4 


IDNINT 


3.5.6 


IFIX 


3.4 


INT 


3.4 


INT2 


3.4 


ISIGN 


3.5.4 


MAXO 


3 . 3 


MAXl 


3.3 


MINO 


3.3 


MINI 


3.3 


MOD 


3.5.1 


NINT 


3.5.6 


REAL 


3.4 


SIGN 


3.5.4 


SIN 


3.2.3 


SINK 


3.2.6 


SNGL 


3.4 


SQRT 


3.2.1 


TAN 


3.2.3 


TANK 


3.2.6 
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APPENDIX 3 

ARITHMETIC TABLES 

TABLE OF POWERS OF T\JO 



2" 


n 


2-n 




1 





1.0 




2 


1 


0.5 




4 


2 


0.25 




H 


:! 


0.125 




IC 


4 


0.0G2 


5 


32 


5 


0.031 


25 


G4 


G 


0.015 


G25 


12H 


7 


0.007 


812 



25G S 0.003 DOG 25 

512 9 0.001 953 125 

1 024 10 0.000 97G 5G2 5 

2 04H 11 0.000 4H.S 2.S1 25 

4 09G 12 0.000 244 140 G25 

« 192 13 0.000 122 070 312 5 

16 3«4 14 O.OOO OGl 035 15G 25 

32 7G.S 15 0.000 03(1 517 57s 1:^5 

G5 53G IG 0.000 015 25S 7X9 0G2 5 

131 072 17 0.000 007 029 394 531 25 

2G2 144 IS 0.000 003 .S14 G97 2G5 G25 

524 2HH 19 O.OOO 001 907 34s G32 .S12 5 

1 04S 57G 20 0.000 000 953 G74 31G 40(1 25 

2 097 152 21 0.000 000 47G .s37 15H 203 125 

4 194 304 22 0.000 000 23S 41S 579 101 5G2 5 

H 3HH GOS 23 0.000 000 119 209 2s9 550 7S1 25 

IC 777 210 24 0.000 000 050 004 044 775 390 025 

33 554 432 25 0.000 000 029 s02 322 3.S7 095 312 5 

G7 lOS SG4 20 0.000 000 014 901 101 193 ^47 05G 25 

134 217 72K 27 0.000 000 007 450 5sO 590 92:) s2S 125 

20S 435 450 2s O.OOO 000 003 725 290 29s 401 914 002 5 

53G M70 912 29 0.000 000 001 S02 045 149 230 957 031 25 

1 073 741 S24 30 0.000 000 000 931 322 574 015 47s 515 025 

2 147 4H3 G4S 31 0.000 000 000 405 GGl 2s7 ;!07 739 257 S12 5 

4 294 9G7 290 32 0.000 000 000 232 s30 (i43 053 s(;9 02s 900 25 

5 5S9 934 592 33 0.000 000 000 110 415 :!21 S2(l 934 sl4 45:) 125 

17 179 SG9 1S4 34 0.000 000 000 OSS 207 000 913 407 407 220 502 5 

34 359 73.S 3GS 35 0.000 000 000 029 103 S30 45G 733 703 013 2.sl 25 

iO 0.000 000 000 014 551 915 22s ;j(iC) s51 sOG ()40 025 

0.00(» 000 000 007 275 957 014 Is;! 425 903 ;',20 312 5 

0.000 000 000 003 G37 97s s07 091 712 051 GOO 15G 25 



GS 


719 


470 


730 


30 


137 


43H 


953 


472 


37 


274 


H77 


900 


944 


3S 


549 


755 


SI 3 


sss 


39 


099 


511 


C27 


770 


40 



0.000 000 000 001 sis 9S9 403 545 s50 475 S30 07s 125 
770 40 0.000 000 000 000 909 494 701 772 92S 237 915 039 002 5 
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HEXADECIMAL ADDITIOI:! TABLE 





1 


2 


3 


4 


5 





7 


s 


9 


A 


B 


C 


D 


E 


F 




i 


2 


3 


4 


5 


li 


7 


N 


9 


A 


B 


C 


D 


E 


F 


10 


1 


2 


3 


4 


5 


6 


7 


s 


9 


A 


B 


C 


D 


E 


F 


10 


11 


2 


"' 


4 


5 


G 


7 


» 


9 


A 


B 


C 


D 


E 


F 


10 


11 


12 


3 


4 


5 


G 


7 


n 


9 


A 


B 


C 


D 


!•; 


F 


10 


11 


12 


13 


4 


5 





7 


H 


9 


A 


B 


C 


I) 


E 


F 


10 


11 


12 


13 


14 


5 


() 


7 


^s 


9 


A 


B 


C 


D 


E 


F 


10 


11 


12 


13 


14 


15 
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Absolute Value Function, 1-2, 3-19 

Accuracy of Function Subprograms, 2-1 
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Argument, Complex, 3-3, 3-5, 3-7, 3-24 

Arithmetic, Floating-Point, 2-1 

Complex Argument, 3-3, 3-5, 3-7, 3-24 
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Complex Functions, 2-2, 3-4, 3-7 

Complex Number, 3-1, 3-3, 3-5, 3-19, 3-20, 3-24 

Complex Numbers, Manipulation of, 1-2, 3-19 

Complex Routines, 2-3/2-4 

Computation, Floating-Point, 2-1 

Conversions, Type, 2-2, 3-16, 3-21 

Description of Function Subprograms, 3-1 
Double-Precision Product of Single-Precision 
Factors Function, 1-2, 3-21 

Elementary Functions 

Exponential Function, 3-3, 3-4 

Hyperbolic Function, 3-11, 3-12, 3-13 

Inverse Trigonometric Function, 3-8, 3-9, 3-10, 3-11 
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Trigonometric Function, 3-4, 3-5, 3-6 
Error 
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Relative Error, 2-1, 3-2, 3-5, 3-7, 3-9, 3-11, 3-20, 3-21 
Error Messages, Generation of, 2-3/2-4 
Error Propagation, 2-1 
Exponential, 1-1, 2-2 
Exponential Function, 3-3, 3-4 
External Versions of Functions, 2-3/2-4 

Floating-Point Arithmetic, 2-1 
Floating-Point Computations, 2-1 
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Index-2 ^8-025 FOO ROO 



INDEX (Continued) 

Positive Difference Function, 1-2, 3-20 

Raising to Powers (Involution), 3-22 
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